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Asymptotic predictions using short-time data in oscillating billiards 
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Faculty of Mathematics and Computer Science, Weizmann Institute of Science, FOB 26, Rehovot 76100, Israe¥] 

Particle motion in a smoothly oscillating non-integrable billiard is known to result in unbounded 
energy growth. Though the asymptotic energy growth rate of an ensemble of particles in an os- 
cillating chaotic billiard is known to be quadratic, there are no estimates available for smoothly 
oscillating pseudo-integrable billiards. The energy growth rate in such systems is so slow that it is 
very hard to predict the asymptotic rates from finite time computations. In this paper, a method 
is proposed to estimate the asymptotic energy growth rate in a system by using short-time data. 
The idea is applied to the case of an oscillating pseudo-integrable system, and it is shown that the 
asymptotic energy growth rate in such systems could be exponential. 



I. INTRODUCTION 



f— *^ For particle motion governed by an autonomous Hamiltonian, the total energy remains constant. If we now add 

^S| a small time-dependent smooth perturbation, the total energy is no longer independent of time. For such systems, 

the energy can remain bounded or grow unboundedly with time. For Id autonomous Hamiltonians, one can show 

using the KAM theorem [TJ [5] that the energy will remain bounded under the influence of a small time-dependent 

p^ perturbation. Though KAM theorem is applicable to 2d and higher dimensional integrable systems too, the energy 

I growth in such systems can still be unbounded due to a phenomenon known as Arnold diffusion [3]. However, this 

H process is exponentially slow and hard to observe in practical situations. For the case of chaotic Hamiltonians, it was 

C/) recently shown that the energy can grow unboundedly, at a much faster rate compared to Arnold diffusion, if a small 

^ time-dependent perturbation is added to the system [4 . 

O I One could carry out a similar analysis for an important class of non-smooth Hamiltonians, namely, the dynamical 

^ billiard. A billiard is a dynamical system in which a particle moves in straight lines within a region bounded by rigid 

O boundaries, and undergoes specular reflections on collision with the boundary jSHZl- Particle motion in oscillating 

' 0/^ billiards has been studied in the context of Fermi acceleration [SHlOj. The concept of Fermi acceleration has found 

|>~» immense applications in areas like coUisional heating in plasma RF sheaths [11] and models of nuclear fission ^12j . 

T^ There also have been several studies of the Fermi accelerator from a quantum mechanical viewpoint |13l [2] . For the 

, ^~ < Id Fermi-Ulam model, it has been proved, using the KAM theorem, that the particle energy remains bounded [15J. 

For the 2d case, it has been numerically shown that energy can grow unboundedly in a smoothly oscillating billiard 

T— I if the frozen system is chaotic jTBH^U]. A theory has also been proposed to explain this observation |21| . 

^ In between completely integrable and completely chaotic systems, there is an interesting class of pseudo-integrable 

Cn systems. Such systems are non-integrable but are not chaotic (all Lyapunov exponents are zero). Liouville- Arnold 

theorem |22| states that if a Hamiltonian system with n degrees of freedom has n independent integrals in involution, 

then it can be integrated by quadratures. However, this notion of integrability is only local. There are systems 

which have the required number of independent integrals in involution, but do not have global action-angle variables 

"n [^ ^^ . And this is because the manifold corresponding to the particle motion in these systems is not diffeomorphic 

to a sphere or a torus, but to a surface of genus greater than one. Rational polygons belong to this interesting class 

of pseudo-integrable systems |25fB7| . A natural question arises regarding the energy growth of a particle if a small 

time-dependent smooth perturbation is added to such systems. Recently, it has been numerically shown that a particle 

experiences unbounded energy growth in smoothly oscillating versions of pseudo-integrable billiards [28j • 

For smoothly oscillating chaotic billiards, it has been numerically observed that the rate of energy growth is 
^ quadratic in time [17J. And this observation is also supported by theoretical arguments [5T1. However, for the case of 

Cu simply connected pseudo-integrable systems, there are no known estimates for the rate of energy growth. This paper 
is the first attempt to give some numerical estimate for the asymptotic energy growth rate in such systems. 

Carrying out numerical computations for energy growth in oscillating billiards is computationally intensive. And 
the lack of accuracy after a certain number of collisions puts a limit on the total time for which the simulation can 
be carried out. This can have important implications for estimating the asymptotic energy growth rate of particles 
in a given system. Basically, if the growth rate up to a certain time has been observed to be quadratic, there is no 
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guarantee that a higher growth rate will not be observed if the computations could be carried out for a longer time. 
Thus, it is very important to have some indirect way of predicting the asymptotic growth rate based on short-time 
data. This paper aims at proposing one such method. 

II. OSCILLATING BILLIARDS 

Consider a 2d billiard with an oscillating boundary. The phase space variables are {x„, y„;M„, w„; t„} where n 
denotes the number of collisions with the oscillating boundary. The particle energy is En — u^ + v'^. In general, u„ 
and v„ are not constant. For the sake of simplicity, let us consider a billiard whose oscillating boundary is just a 
straight line and moves in a direction normal to itself. Let the wall velocity at the instant of the nth collision be Wn 
and always along the y-axis. The change of velocity of the particle in the nth collision with the oscillating boundary 
is 

Au„ = 2w„ (1) 

The time interval between two collisions with the oscillating boundary is 

Ai„ - -^ (2) 



where Z„ is the distance travelled by the particle in between the nth and (n — l)th collision. As particle energy grows 
unboundedly {E^ — > oo), the time-interval between collisions decays to zero (Ai„ — )■ 0), and in this limit, we have 

^" = Ve'-^ (3) 



dt I 



Substituting v — \fE cos (/), we get 
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(4) 



Now, what we are interested in is the growth rate averaged over the time period of oscillations of the base. Time 
averaging Eq. Q, we get 

dE — / d(b 'iw \ 

dt \ dt I cos (p / rp 

where T — 2tt/uj is the time period of oscillations of the billiard. For the Id Fermi-Ulam model, cos 4> = I, and I is a 
known periodic function if E is high enough. Thus, in this case, the average of the right hand side of Eq. ([5]) is zero 
and there is no net acceleration of the particle. However, for a general 2d billiard, I and cos cf) change randomly, and 

if the average, ( 2 tancj)-^ + , "^^ , )> 0, then this leads to an unbounded growth of energy. If the long time average, 
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dt^R>0 (6) 



then the particle also experiences exponentially fast growth of energy, E — Eqc^'* . However, for a given system, if R 
is very low, then it is hard to see the asymptotic exponential growth numerically. This is mainly because numerically 
accurate computations can be done only for a finite time, which may be too short to observe the exponential energy 
growth in some systems. In this paper, a way of estimating asymptotic growth rates using short-time data is proposed. 

III. MAKING ASYMPTOTIC PREDICTIONS 

The most general equation for evolution of energy of an ensemble of particles in any system is 

f=/(i=,0 (7) 



where / is an arbitrary function of E,t. The method proposed in this paper is appHcable to a sub-class of Eq. (|7|, 
namely, those equations where / {E,t) is variable separable, 



dE 
~dt 



f{E,t)^g{E)h{t) 



(8) 



where g{) and h{) are two arbitrary functions of their arguments. As can be seen from Eq. (pi), the oscillating 
billiard certainly belongs to this class of systems. If both the functions g, h are known analytically, then one can 
surely go ahead and do rigorous analysis. However, in certain situations like particle motion in oscillating billiards, 
the function, hit), has terms that change in a way that is random for all practical purposes. In such situations, one 
cannot do any rigorous analysis for predicting the system behavior. Also, as mentioned earlier, one also cannot do 
computations for very long times due to limits on the accuracy. 

The basic idea proposed in this paper is that measuring the dependence of short-time growth rates on the initial 
energy of the system can give us some indication of asymptotic behavior. For various asymptotic energy growth rates, 
it can be shown, using Eq. dsl, that the dependence of the linear part of the growth rate on initial energy, _Eo, must 
be 



Linear growth rate 

Quadratic growth rate 

Exponential growth rate 



E = Eo + at 

E = Eo + f3E°-h 

E = Eo + -fEot + 



(9) 



where a, (3,^ are arbitrary constants. Thus, the idea is to carry out the simulations till a time for which the growth 
rate is predominantly linear and then, using Eq. (pi, we can say that 

1. If a plot of {E — Eq) is independent of Eq, the asymptotic growth rate is linear 

2. If a plot of {E — Eq) / Eq-^ is independent of Eq, the asymptotic growth rate is quadratic 

3. If a plot of {E — Eq) I Eq is independent of Eq, the asymptotic growth rate is exponential 

Now, we use this idea to estimate the asymptotic energy growth of an ensemble of particles in an oscillating trapezium 
(see Fig. ITl, which is a pseudo-integrable system. For the case of a trapezium with an oscillating base. Fig. p] shows 
that {E — Eq) j Eq is independent of Eq (approximately), hence indicating that the asymptotic energy growth rate in 
this billiard could be exponential. In Fig. pi a plot of (E — Eq) /Eq-^ can be clearly seen to depend on Eq. Though 
the plots in Figs. |2]and|3]are not conclusive, they do provide a reasonable indication of the asymptotic energy growth 
rate being exponential. 




Figure 1: Geometry of the trapezium. 
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Figure 2: Plot of {E — Eo) /Eq vs. time for an ensemble of 1000 particles in an oscillating trapezium with parameters: L — 0.5, 
h = 4.0 and ZADC = O.Stt + 7r/l8. Amplitude of oscillations of the base (AB) was 0.1 with a frequency of lj = 0.1, and with 
an initial horizontal velocity, wo = ^ui/n . 
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Figure 3: Plot of {E 
as in Fig. [2] 



Eo) /Eq'^ vs. time for an ensemble of 1000 particles in an oscillating trapezium with the same parameters 



IV. CONCLUSIONS 

Though the idea presented in this paper has been appHed to the specific case of an oscillating biUiard, it is very 
general in nature. There are many problems in engineering and applied sciences where prediction of asymptotic 
behavior is of crucial importance. One limitation of the method proposed in this paper is that, using this method, 
it is hard to differentiate between two growth rates that are not well separated, eg. exponential versus a polynomial 
of high order. However, I hope that the idea proposed in this paper will provide a starting point in developing more 
elaborate tools to handle this problem. 

Acknowledgments 

I acknowledge support of the Israel Science Foundation and the Minerva Foundation. I would also like to thank 
Vered Rom-Kedar and Dmitry Turaev for interesting discussions. 



[1 
[2 
[3 
[4 
[5 
[6 

[7; 

[8 
[9 

[10 

[11 
[12 
[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 
[22 
[23 
[24 
[25 
[26 
[27 
[28 



M. Tabor, Chaos and Integrability in Nonlinear Dynamics: An Introduction (New York, Wiley, 1989) 

M. A. Lieberman and A. J. Lichtenberg, Phys. Rev. A 5, 1852 (1972) 

V. I. Arnold, Sov. Math. Doklady 5, 581 (1964) 

V. Gelfreich and D. Turaev, Conim. Math. Phys. 283, 769 (2008) 

G. D. Birkhoff, Acta Math. 50, 359 (1927) 

Ya. G. Sinai, Russian Math. Surveys 25, 137 (1970) 

L. A. Bunimovich, Commun. Math. Phys. 65, 295 (1979) 

E. Fermi, Phys. Rev. 75, 1169 (1949) 

S. M. Ulani, On some statistical properties of dynamical systems. Proc. 4th Berkeley Symp. Mathematical Statistics 

Probability Univ. Cahfornia 3, 315 (1961) 

G. M. Zaslavsky, Chaos 15, 015103 (2005) 

M. A. Lieberman and V. A. Godyak, IEEE Trans. Plasma Sci. 26, 955 (1998) 

J. Blocki, Y. Boneh, J.Nix, J.Randrup, M.Robel, A.Sierk and W.Swiatecki, Annals Phys. 113, 330 (1978) 

A. J. Makowski and S. T. Dembinski, Phys. Lett. A 154, 217 (1991) 

J. V. Jose and R. Cordery Phys. Rev. Lett. 56, 290 (1986) 

L. D. Pustyl'nikov, Theo. Math. Phys. 57, 1035 (1983) 

J. Koiller, R. Markarian, S. O. Kamphorst and S. P. de Carvalho, Nonlinearity 8, 983(1995) 

A. Loskutov, A. B. Ryabov and L. G. Akinshin, J. Phys. A: Math. Gen. 33, 7973 (2000) 

E. D. Leonel and P. V. E. McClintock, J. Phys. A: Math. Gen. 38, 823 (2005) 

R. E. de Carvalho, F. C. de Souza and E. D. Leonel, J. Phys. A: Math. Gen. 39, 3561 (2006) 

S. O. Kamphorst, E. D. Leonel and J. K. L. da Silva, J. Phys. A: Math. Theor. 40, F887 (2007) 

V. Gelfreich and D. Turaev, J. Phys. A: Math. Theor. 41, 212003 (2008) 

V. I. Arnold, Mathematical Methods of Classical Mechanics (New York, Springer, 1989) 

V. V. Kozlov, Russian Math. Surveys 38, 1 (1983) 

J. J. Duistermaat, Comm. Pure. Appl. Math. 33, 687 (1980) 

P. J. Richens and M. V. Berry Physica D 2, 495 (1981) 

E. Gutkin, Physica D 19, 311 (1986) 

E. Gutkin, J. Stat. Phys. 83, 7 (1996) 

K. Shah, D. Turaev and V. Rom-Kedar, Phys. Rev. E 81, 056205 (2010) 



